Stochastic dynamics of electronic wave packets in fluctuating laser fields 
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The dynamics of a laser-excited Rydberg electron under the influence of a fluctuating laser 
field are investigated. Rate equations are developed which describe these dynamics in the limit of 
large laser bandwidths for arbitrary types of laser fluctuations. These equations apply whenever 
all coherent effects have already been damped out. The range of validity of these rate equations is 
investigated in detail for the case of phase fluctuations. The resulting asymptotic power laws are 
investigated which characterize the long time dynamics of the laser-excited Rydberg electron and it 
is shown to which extent these power laws depend on details of the laser spectrum. 
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I. INTRODUCTION 



The advancement of sophisticated trapping techniques 
and the development of powerful laser sources has stimu- 
lated numerous theoretical and experimental investiga- 
tions on the dynamics of wave packets in elementary 
quantum systems. An understanding of their dynamics 
is important for our conception of quantum mechanics 
and its relation to classical mechanics. So far most of the 
research in this context has concentrated on coherent as- 
pects of wave packet dynamics which may be traced back 
theoretically to semiclassical aspects originating from the 
smallness of the de Broglie wave lengths involved. How- 
ever, for an understanding of the emergence of classical 
behaviour also a detailed understanding of the destruc- 
tion of quantum coherence is required. Such a destruc- 
tion of coherence may arise from the coupling of a quan- 
tum system to a reservoir or from stochastic external 
influences. Though by now many aspects of the coherent 
dynamics of wave packets are well understood still many 
questions concerning the influence of stochastic pertur- 
bations on elementary quantum systems with a high level 
density are open. 

A paradigm of a quantum system in which many of 
these latter aspects can be investigated in great detail 
are Rydberg systems interacting with fluctuating laser 
fields. Due to the inherent stochastic nature of laser light 
a detailed understanding of optical excitation processes 
with fluctuating laser fields is of vital interest for laser 
spectroscopy. Rydberg systems jl|-[| are of particular 
interest in this context due to their high level density 
close to an ionization threshold. Any laser-induced ex- 
citation process which involves Rydberg and continuum 
states close to an ionization threshold typically leads to 
the preparation of a spatially localized electronic Ryd- 
berg wave packet |Q . Under the influence of a fluctuat- 
ing laser field the coherence of such an electronic wave 
packet is disturbed whenever it is close to the ionic core 
where the electron- laser interaction is localized || • Even- 
tually these random perturbations are expected to lead 



to a stochastically dominated Brownian motion of the 
electronic wave packet. Indeed it has been demonstrated 
recently || that such a transition to a diffusive behaviour 
takes place and that this diffusive dynamics are domi- 
nated by characteristic power laws which govern the time 
evolution of the Rydberg system. However, these previ- 
ous investigations were restricted to a particular type of 
phase fluctuations of laser fields which can be described 
by the so called phase diffusion model (PDM) This 
PDM implies a Lorentzian spectrum for the fluctuating 
laser field. It is known from the dynamics of three level 
systems that the somewhat unrealistic asymptotic fre- 
quency dependence of a Lorentzian laser spectrum may 
lead to unphysical predictions ||. However, from our 
previous investigations it remained open to which extent 
these characteristic, diffusive long time dynamics of an 
excited Rydberg electron depend on details of the fluc- 
tuations of the exciting laser field. Such a dependence 
might be expected on intuitive grounds as the diffusion 
of the excited Rydberg electron in energy space eventu- 
ally also reaches the far-off resonant regions of the laser 
spectrum. 

In this paper we tackle these open questions by gen- 
eralizing our previous results to arbitrary types of laser 
fluctuations. For this purpose we derive rate equations 
for the relevant density matrix elements of the excited 
Rydberg electron which are averaged over the fluctua- 
tions of the laser field. These rate equations are based 
on a decorrelation of the relevant electron-field averages. 
This decorrelation approximation (DCA) is valid as long 
as the characteristic correlation time of the fluctuating 
laser field, i.e. its inverse bandwidth, is much smaller 
than all other relevant intrinsic dynamical time scales. 
Within this framework it will become apparent that it is 
the laser spectrum only which determines the time evo- 
lution of the excited Rydberg electron. On the basis of 
this approach it will be demonstrated which aspects of 
the diffusive long time dynamics of an excited electronic 
Rydberg wave packet depend on which details of the laser 
spectrum. 

The range of validity of these Pauli-type rate equa- 
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tions is investigated in detail for a special class of phase 
fluctuations P,|10| of the exciting laser field. This special 
class of phase fluctuations implies non-Lorentzian spectra 
which for large laser frequencies decrease more rapidly 
than a Lorentzian. These phase fluctuations might be 
considered as a realistic model for a single mode laser 
field which is operated well above the laser threshold. In 
order to access the range of validity of the Pauli-type 
rate equations for this special class of laser fluctuations a 
more general master equation is derived for the averaged 
density operator of the Rydberg electron. This more gen- 
eral approach is also capable of dealing with all coherent 
aspects of the laser excitation process. 

The paper is organized in the following way: In Sec. 
II the theoretical models for describing the laser fluc- 
tuations and the dynamics of the Rydberg electron are 
presented. For the sake of simplicity we restrict our sub- 
sequent discussion to Rydberg systems which can be de- 
scribed within the framework of a one-channel approxi- 
mation Typically this approximation is well satis- 
fied for Alkali atoms. In Sec. Ill we derive rate equations 
which describe the dynamics of the excited Rydberg elec- 
tron averaged over the laser fluctuations. Self consistent 
validity conditions for the applicability of the decorrela- 
tion approximation (DCA) are discussed on which these 
rate equations are based on. In Sec. IV characteristic 
aspects of the time evolution predicted by the rate equa- 
tions of Sec. Ill are exemplified. The different dynamical 
long time regimes and their characteristic power laws are 
discussed in detail. From the resulting analytical expres- 
sions for these power laws it is apparent to which extent 
they depend on details of the laser spectrum. In Sec.V we 
derive a more sophisticated master equation for the av- 
eraged dynamics of the excited Rydberg electron. This 
master equation is capable of describing also coherent 
aspects of the dynamics of the excited Rydberg electron 
but its validity is restricted to a particular class of phase 
fluctuations only. In Sec. VI solutions of this master 
equation are compared with the corresponding results of 
the rate equations. Thus we are able to determine the 
range of validity of the DCA. In our subsequent discus- 
sions we use atomic units (m e = e = U = 1). 



II. THEORETICAL FRAMEWORK 

In this section the theoretical models are introduced 
with which the fluctuating laser field and the excited Ry- 
dberg system are described. 



We assume that this laser field can be described by a 
classical stochastic process [0. The mean frequency of 
this laser field is denoted u> and e is its polarization vec- 
tor. The fluctuations of this laser field are described by 
the envelope function e(t) which is assumed to be slowly 
varying on time scales of the order of 1/uj. The associated 
spectrum of this laser field is defined by H 
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S(Q) = -Re / dTK{r)e- inT 
n Jo 



(2) 



with the two-time correlation function of the slowly vary- 
ing amplitude 



K(T) = (e{t + r)e*(t)) 



(3) 



Thereby (...) denotes statistical averaging over the fluc- 
tuations of the laser field. 

For a single mode laser which is operated well above 
laser threshold to a good degree of approximation the 
amplitude of the laser field is stable. In these cases fluc- 
tuations of a realistic laser field can be described by a 
classical electromagnetic field whose phase is fluctuating, 
i.e. 



e(t) = e e~ 



Ht) 



(4) 



The fluctuating phase $(t) obeys the (Ito-)stochastic dif- 
ferential equation 



(5) 



d<t>{t) = -4>(t)/3dt + V2b(3dW(t) 



with cj)(t) = &(t). In Eq.(||) determines the correla- 
tion time of the stochastic frequency <j>(t) and dW(t) is 
the differential of a real- valued Wiener process with zero 
mean and unit variance, i.e. (dW(t)) = 0, (dW(t) 2 } = dt 
JlT|. Eqs.ffl) and (§) imply the relation 



K(t) = \e \ 2 exp [-6t + b/(3 (l - e~ r/3 ) 



(6) 



In the limit of large values of j3 $(£) itself approaches a 
real- valued Wiener process, i.e. 



(7) 



d<f>(t) = V2bdW{t) 



This limiting case constitutes the so called phase diffusion 
model (PDM). It implies a Lorentzian laser spectrum of 
the form 



S(Q) 



M 2 - 



7T b 2 + n 2 ' 



(8) 



A. The fluctuating laser field 

We consider an atomic or molecular Rydberg system 
which is driven by a laser field with electric field strength 



E(t) = ee(t)e 



c.c. 



(1) 



For (3 ^> b the parameter /3 may be interpreted as a cut- 
off parameter of the laser spectrum. This becomes appar- 
ent by noting that for frequencies f2 <C j3 the spectrum is 
always approximately Lorentzian whereas for large fre- 
quencies, i.e. f2 3> j3, it tends to zero more rapidly as 
can be seen in Fig.ll]. More precisely, for /3 3> b we obtain 
from Eqs.ra) and (0) the approximate relation 
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the initial state \g) and an energy normalized continuum 
state |e) with energy e w 0. Within the framework of 
QDT a as well as d eg are approximately energy indepen- 
dent close to the ionization threshold, i.e. for energies 
|e| <C 1. A one channel approximation is appropriate for 
all cases in which excited states of the ionic core are lo- 
cated far away from the excited energy region. Typically 
this condition is fulfilled for Alkali atoms. 

The main problem is to solve the stochastic 
Schrodinger equation associated with the Hamiltonian of 
Eq. ( |l0| ) . In general this is a complicated task due to the 
simultaneous presence of the laser fluctuations and of the 
threshold effects arising from the infinitely many bound 
Rydberg states converging to the ionization threshold. 
By the resulting intricate interplay between laser fluctu- 
ations and threshold phenomena it is difficult to apply 
stochastic simulation methods which typically become 
unreliable due to numerical inaccuracies in particular for 
long interaction times. 



FIG. 1. Lorentzian laser spectrum according to Eq.(g) 
(solid line) and non-Lorentzian spectrum with f3 = 206 
according to Eqs.(^) and (|^) (dashed line). 



B. The interaction Hamiltonian 

Let us assume that the considered atom or molecule 
is prepared initially in an energetically low lying bound 
state |<?) with energy e g . Furthermore, it is excited res- 
onantly by the fluctuating laser field to Rydberg and/or 
continuum states \n) with energies e n close to an ioniza- 
tion threshold. In the dipole and rotating wave approx- 
imation this excitation process can be described by the 
Hamiltonian 



H= £ |iXi|ei-E[ £ (*) e ^ 



h.c] . (10) 



]=B,n 



The dipole matrix elements between the initial state \g) 
and the excited Rydberg states \n) are denoted d ng — 
(n\d • e \g). It is understood that the sum over the ex- 
cited states \n) appearing in Eq.(|i"o|) includes also an inte- 
gration over the adjacent continuum states. The energies 
of the excited Rydberg states and the energy dependence 
of the dipole matrix elements entering the Hamiltonian 
of Eq.(hd) can be determined with the help of quantum 
defect theory (QDT) jij. In the case of excited Rydberg 
states which can be described within the framework of 
the one-channel approximation we find, for example, 
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2(n-a) 2 ' 



d ng = d eg (n - a) 3/2 = d eg \t 



1/2 



(11) 



Thereby a denotes the quantum defect of the excited Ry- 
dberg series and d eg is the dipole matrix element between 



III. DECORRELATION APPROXIMATION 
(DCA) 

In this section an approximation method is developed 
for determining the dynamics of the Rydberg system in 
the fluctuating laser field. This approximation method 
is based on a decorrelation of atom-field averages and 
leads to a Pauli-type master equation for the density op- 
erator of the Rydberg system which is averaged over the 
laser fluctuations. In this master equation all coherences 
between different energy levels have been eliminated adi- 
abatically. This decorrelation approximation (DCA) is 
valid for arbitrary types of laser fluctuations provided 
these fluctuations arc sufficiently fast (compare with con- 
ditions ©, |L|) and (|l|) derived below). 

Let us start by determining first of all an approxi- 
mate equation of motion for the probabilities p nn (t) = 
(n\ip(t))(ip(t)\n) of observing the excited Rydberg sys- 
tem in one of the Rydberg states \n). Neglecting coher- 
ences p n n'(t) with n 7^ n 1 we find from the stochastic 
Schrodinger equation with Hamiltonian ( |l0| ) the relation 

Pun (*) = 2 | d ng | 2 Re{ f dt'e^e*^-^-^-^} 



x [Pgg(t') -Pnn{t')] 

- 'id ng e(t)e^^ (t ^p gn (to) + h.c. (12) 



with < tn 



< 



t 



and with Pkj(t) = 
(k\ip{t)) {ip{t)\j) exp[ii(efc — ej)]. The mean excited energy 
is denoted e = e g + uj + Soj with the quadratic Stark shift 
contribution of all other (non-resonant) states. In the 
subsequent discussion we assume for the sake of simplic- 
ity that the intensity dependence of 8ui does not affect 
the dynamics of the excited Rydberg system. This is 
valid either for pure phase fluctuations of the laser field 
or in the case of arbitrary laser fluctuations for laser 
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bandwidths which are much larger than Slo. The self- 
consistency condition for the omission of the coherences 
Pnn' (t) with n =/= n' will be discussed later (compare 
with Eq.(|l9|)). Taking the integration interval [to,t] to 
be smaller than the characteristic time scale over which 
[Pnn{t') — Pgg(t')] changes significantly this latter term 
can be approximated by its value at time t. If on the 
other hand the interval [to,t] is assumed to be larger than 
the correlation time of the fluctuating laser field we can 
replace the lower integration limit to by — oo in Eq.(12). 
Thus we obtain 



P nn (t) = 2\d ng \ 2 \figg(t) ~ p n n{t)\ 
poo 

xRe / dTs(t)e*(t- t} 



(13) 



,-i-r(e-e„) 



-[idngE^e-^- 1 ^-^ Pgn (t ) + h.C.}. 



Now we are able to carry out the statistic average (...) 
over the laser fluctuations. Due to the above mentioned 
conditions on the integration interval the involved density 
matrix elements and the laser field e(t) decorrelate. As 
(e(t)) = the contribution of the last term on the right 
hand side of Eq. (|l3[) vanishes. The remaining terms yield 
the rate equation 



(Pnn {t)) =n ng [( Pgg (t)) - (pnn(*)> 

with the time independent rates 

K ng = 2ir\d n g\ 2 S(e - e n ) 



(14) 



(15) 



and with the laser spectrum S(Sl) as defined by Eq.(g). 

In order to work out quantitative criteria for the valid- 
ity of Eq.(|l4|) let us define an effective bandwidth B of 
the fluctuating laser field S(tt) by the relation 



BS(0)tt 



dns(n) = <|e| 2 ) 



(16) 



The quantity 1/B measures the correlation time of the 
fluctuating laser field. In the case of the PDM, for ex- 
ample, this effective bandwidth equals the width of the 
Lorentzian spectrum b. In all other cases which are de- 
scribed by Eq.(^|) it characterizes the effective frequency 
width of the laser spectrum of Eq.(||). According to 
Eq.(|l4|) the inverse rates l/lZ ng define the characteris- 
tic time scale over which [p gg (t) — p n nit)] varies signifi- 
cantly. Thus, the above decorrelation approximation ap- 
plies to cases only for which B ^> lZ ng - Typically one 
finds lZn rc3 g > TZ-ng with e„ rcs = e so that one of the va- 
lidity conditions for the decorrelation condition becomes 



B > 27r|d„ 



2 5(0). 



(17) 



Eliminating S(0) by Eq.(|l6|) we finally arrive at the 
equivalent condition 
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b 2 » - (n R 



(18) 



with the average Rabi-frequency (Qr) = 
2|c?„ rcsff |-\/(|e| 2 ). What remains to be found is a valid- 
ity condition for neglecting the coherences p nn '{t) with 
n =/= n' in Eq.(|l2|). As long as 



s +l 



{Hi 



2B 



(19) 



the coherences p nn ' (t) with n ^ n' are rapidly oscillating 
functions in comparison with the slowly varying probabil- 
ities (pnn(t)) and (p gg (t)) entering Eq. (|l2|) so that their 
influence averages to zero approximately. Therefore the 
inequalities dl8| ) and ( |l9| ) are the required conditions for 
the validity of the DCA. According to these conditions 
we may distinguish two limiting cases. In the limit of 
small laser bandwidths for which B <^ \(-n rc3 — en,.,, s +i| 
they reduce to the requirement B 2 3> (Or) /2. In two- 
level systems which are excited resonantly by a fluctuat- 
ing laser field this is the well known limit of large laser 
bandwidths in which the dynamics are dominated by rate 
equations. In the opposite limit where the bandwidth is 
large enough to affect many excited Rydberg states, i.e. 
B ^S> \t-n rcs — £ri rcs +i|, the conditions for the applicability 
of the DCA reduce to the relation B ^ 7/t. Thereby we 
have introduced the laser-induced rate 



7 = 2^K 9 | 2 (M 2 )^ - (Q R ) 2 \e r , 



e^+xl" 1 (20) 



which characterizes ionization of the initial state \g) into 
continuum states close to the ionization threshold accord- 
ing to Fermi's Golden rule. 

Up to now, our arguments for the derivation of the 
rate equation (f|) and of conditions © and (£f) apply 
for discrete excited states only. However, our previous 
arguments can be generalized easily also to continuum 
states by viewing these continuum states as infinitesi- 
mally spaced discrete energy levels. According to quan- 
tum defect theory close to an ionization threshold the 
energy dependence of the discrete dipole matrix elements 
d ng is described by Eq. (|ll|) . Thus condition (|l^) reduces 
to 



B»m a \ 2 (\e\ 2 )^-. 



(21) 



In the limit of an infinitcsimally small level spacing be- 
tween the excited states Eqs. (|l4j) and ([ll]) imply that 
the probability of finding the excited Rydberg system in 
a continuum state becomes vanishingly small. Thus we 
find 



{Pee{t))=n eg ( Pgg (t)) 



(22) 



with K eg = 2ir\d eg \ 2 



S(e — e). Integration over the whole 
electron continuum finally yields 



(Pton (t))=T( Pgg (t)) 



(23) 
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with the mean ionization probability (Pj orl (i)) = with 
L°° de(p ee (t)} and with the effective ionization rate 

/>oo 

r = 2vr / de\d eg \ 2 S{e- e). (24) 
Jo 



If the mean excited energy e is located well above thresh- 
old and if d eg is still energy independent over the energy 
region over which S(e — e) is significant, this effective ion- 
ization r reduces to the previously introduced ionization 
rate 7 of Eq.©. 

The rate Eqs. ( p"4| ) and (|2^) together with the conser- 
vation of probability, i.e. 



(pgg (*)> = 1 - (iW*)) - £>»n(*)> 



(25) 



and together with the initial condition /j(t = 0) = |g)(<7| 
determine the time evolution of a laser excited Rydberg 
electron within the framework of the DCA. 



IV. STOCHASTIC DYNAMICS OF RYDBERG 
SYSTEMS WITHIN THE DCA 

In this section the dynamics of a Rydberg system is 
investigated with the help of the DCA on the basis of 
Eqs.Q, @ and ©. 

Within the framework of a one-channel approximation 
the excited energies and the energy dependence of the 
relevant dipole matrix elements of a Rydberg system can 
be described by Eqs. (pi]). They are characterized by a 
quantum defect a and by an energy-normalized dipole 
matrix element d eg which are both approximately energy 
independent for |e| <C 1. The laser-induced coupling be- 
tween the initial state \g) and the excited Rydberg- and 
continuum states is characterized by the ionization rate 
7 of Eq.(|2(]). Typically this description is adequate for 
Rydberg states of Alkali atoms. 

The rate equations for the averaged density operator 
of the Rydberg system (compare with Eqs.(|l4|),(|23|) and 
(p5|)) can be analyzed in a convenient way with the help 
of Laplace transformations. Defining the Laplace trans- 
formed density operator by 



(P(z)) 



dte lzt (p(t)} 



(26) 



the associated inverse transformation is given by 

n pOQ-\-iO 

(p(t)) = - / dze-^(p(z)). (27) 

Thus the Laplace transformed rate equations ([hi]), ( p^ ) 
and (|25|) imply the relations 



(Pgg( z )) = 
and (P ion (z)) = 



1 



T — iza{z) 
iT 

z [T — iza(z)] 



(28) 
(29) 



(30) 



The rates TZ ng entering Eq.(p0|) characterize the transi- 
tions between states \g) and \n) within the DCA and 
are defined by Eq.(|l5|). In the derivation of Eqs.(|2"g|) 
and ( |29| ) (p gg {t)) has been neglected in comparison with 
EJpnnit)} and (P lon {t)) in Eq.©. 

We may distinguish various dynamics regimes which 
are treated subsequently. 



A. Asymptotic long time behaviour 



The time evolution of the averaged density operator of 
the Rydberg system can be obtained from Eqs.(p8|) and 
( ^9| ) and from the inversion formula (|27|). In general the 
time evolution will exhibit both exponential decays orig- 
inating from poles of the Laplace transforms (£8|) and 
( p9| ) in the complex z-plane and power law decays which 
originate from cut contributions starting from the branch 
point of a(z) at z = 0. As the asymptotic long time be- 
haviour will be dominated by power law decays we have 
to investigate the structure of the characteristic kernel 
a(z) around the branch point z = in more detail. From 
Eqs.© and @ it follows that for z -> its mam con- 
tributions arise from the infinitely many Rydberg states 
very close to the ionization threshold. Hence in the long 
time limit we may approximate S(e — e n ) by S(e) in ex- 
pression ( |l5|) . Furthermore we may replace the sum over 
all Rydberg states in Eq.(|30|) by an integration. So finally 
in the limit z —> we obtain the relation 



a(z) 



2tt fi2Tr\d eg \ 2 S(e) 



3^3 



1/3 



(z-0). (31) 



Inserting Eq.©) into Eqs.©), ©) and ©) one obtains 
the asymptotic long time behaviour 



(Pgg(t)) = 



S(e) 



1/3 r( 5\ 

[ ~ 3) (*y)" 



•5/3 



3(r/ 7 ) 5 



(Pion(t)) = 1 - 



S(e) 



1/3 



r(I) 

3(T/ 7 ) 



-2/3 



(32) 
(33) 



with T(x) = J °° du e~ u u x ~ 1 denoting the gamma func- 
tion [|l3|. Eqs.(|3^) and ( |3^ ) describe the time evolution 
of the mean initial state probability and of the mean ion- 
ization probability for sufficiently long interaction times. 
They are generalizations of our previous results of Rcf. 
H which were only valid for phase fluctuations of the 
PDM. Within the framework of the DCA these asymp- 
totic laws are valid for arbitrary fluctuations of the laser 
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field provided B 3> j/tt (compare with Eq.(|2lJ)). Obvi- 
ously this asymptotic long time behaviour is independent 
of the quantum defect which characterizes the influence 
of the ionic core of the Rydberg system. Furthermore, 
the characteristic exponents of the long time behaviour 
are a peculiar property of the Coulomb problem and do 
not depend on details of the laser spectrum. However, 
the time independent pre-factors of these power laws de- 
pend on S(Q) and on the effective ionization rate T of 
Eqs.(|) and (|J). 

After which interaction time do we expect the asymp- 
totic power laws of Eq s.(|32| ) and ([53]) to become valid? 
As apparent from Eq.(33) this diffusive long time dy- 



namics finally leads to complete ionization of the Ry- 
dberg system. Thus it is reasonable to characterize the 
onset of this asymptotic long time dynamics by a stochas- 
tic ionization time t c which is defined by the condition 
(PioniQ) = 1/2 which yields 



7r(§) 



8 S(e) 

27 w: 



1/2 



(34) 



B. Intermediate interaction times 

In this section we deal with characteristic aspects of 
the dynamics described by Eqs.(|l4|) and ( ^3|) in cases in 
which the interaction times are large enough so that the 
initial state \g) is depleted significantly but which are 
still much smaller than the stochastic ionization time t c 
of Eq.(|34|). For these interaction times we may distin- 
guish two characteristic regimes depending on whether 
the excited states are close to the ionization threshold or 
above or whether they are located well below this thresh- 
old. 



1. Excitation at or above threshold 

In this dynamical regime the significantly excited 
states are located at the ionization threshold or above. 
This implies that all rates lZ ng which describe the cou- 
pling between states \g) and \n) are small in comparison 
with the total ionization rate 7. Thus considering in- 
teraction times t which are not very much larger than 
I/7 implies that we may take tlZ ng <C 1 for all quan- 
tum numbers n. Hence, considering the formal solution 
of Eq.pl) 



(Pnn(t)) = Hr, 



dse 



-(t-s)K ng 



(pgg( s )) 



(35) 



we may replace the exponential by unity. Performing the 
summation over all Rydberg states we find with the help 
of J2 n ^ n 9 ~ 7 — T (the sum over all Rydberg states has 
been replaced by an integral) 



X>nn (t)) = (7 - r) ( Pgg (t)) - 7 2 A Sp / ds (p gg (s)) 

(36) 

where we used the quantity 
A Sp = 



-J-2 / de n (^-)S 2 (e- e„) 
rl 2 ) J-00 dn 

deS 2 (e + e)(2e) 3 / 2 



(37) 



which characterizes the spectral properties of the fluctu- 
ating laser field. 



We take the time-derivative of En. (25) and eliminate 
the ionization probability with Eq. (|23| ) . Inserting Eq. (|3^ ) 
in this equation we find an integro-differential equation 
for (p gg (t)}, namely 

7 (p gg (t)) ~ 7 2 As P / ds ( Pgg (s)) + ( Pgg (<)) = 0. (38) 
Jo 

As As p <C 1 one finally obtains the relations 

(Pm(*)> = 
and 
(Pion{t)) 



1 



1 + A 



Sp 



exp(-Tt) + A Sp exp(7A Sp i) (39) 



7(1 + A Sp ) 



exp(7A Sp i) - exp(-Tt) . (40) 



These equations even apply to interaction times t < I/7. 
Note that consistent with our approximations the inter- 
action times always fulfill the inequality "fAs p t < 1. In 
the special case of the PDM A$ p reduces to 



A W^M( l - 2 wH)- 



(41) 



Very close to threshold, i.e. for e — > 0, one obtains 
Asp = v / ^/(27r) so that in this special case we obtain 
again our previous results of Ref. 



2. Excitation well below threshold 

If the fluctuating laser field excites Rydberg states well 
below the ionization threshold, i.e. B <JC |e| and e < 0, the 
considerations of Sec.lVBl have to be modified. Since 



the interaction time is assumed to be smaller than the 
stochastic ionization time t c we can approximate the 
characteristics of the dominantly excited Rydberg states 
by 



"ng 



e + 2ix{n ~ n res )/Te, 
d eg ^/2tt/T^. 



(42) 
(43) 



G 



Thereby Tj = 2tt(— 2e) -3 / 2 denotes the classical Kepler 
period of the mean excited Rydberg state of energy e < 0. 
Whereas in the previous subsection the stochastic influ- 
ence could be characterized by a single average spectral 
property for arbitrary types of fluctuations, namely by 
As p , the excitation dynamics well below threshold turns 
out to be much more sensitive to the details of the laser 
spectrum. This is easily demonstrated by considering 
phase fluctuations which can be described by the laser 
spectrum of Eq.(||) as a particular example. This spec- 
trum describes fluctuations of a single mode laser field 
well above the laser threshold in the limit j3 ^> b. In 
addition, if only Rydberg states well below threshold are 
excited significantly we may neglect the effective ioniza- 
tion rate T in the denominator of Eqs.(|2§l) and (|29|). In 
the limit 8 3> b we thus arrive at the relations 



(p 99 (t)) = T^r 



(Pionit)) = ~ 

7T 



1 



2 7 6 



A-arctanfA 



2 7 6 



(44) 
(45) 



Thus within this limit for arbitrary values of (3 and b 
the influence of the phase fluctuations of the laser field 
is described by the single scaling function /(r) which is 
defined by the equation 



df(r) 
dr 



= -Im/ ?« 



-iCr 



dx 



1 - i((x 2 + x A ) 



.x 4 ) } ' 



(46) 



To end up with Eq.(|4(|) we had to apply the further 
approximation Q b in the laser spectrum of Eq.(^|). 
Physically speaking this approximation means that we 
consider cases in which the essential dynamics are domi- 
nated by energy states which are located in the wings of 
the laser spectrum. 

In the limits r <C 1 and r ^> 1 asymptotic expressions 
are easily obtained from Eq. ( filf ) . The limit of small val- 
ues of t is realized in the PDM where (3 — > 00 and where 
the spectrum of Eq. (^J) reduces to a Lorentzian form. In 
this case one obtains the expression [0 



f(r) 



(r « 1). 



Consequently Eqs.(|44[) and ( [451) yield 

(Pg a (t)) = 



(Pion(t)) = 2 



■wbTrft 

H— (I 




From the numerical data shown in Fig.|2| it is apparent 
that Eqs.(f48|) and ( |4g| ) are good estimates for interaction 
times t < t PDM = T 7 (3 2 / (200b-f) . 
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FIG. 2. Numerical solutions of Eqs.(|4J) and (|45|) 

a) initial-state probability 

p gg> scaie{t) = (pgg(t)) \T T (3 (solid line) together with the 
asymptotic behaviour according to Eq. mm (dashed) and 
Eq.@ (long dashed). 

b) ionization probability 

Pian,scalS) = *"[/V|e| - arctanCfl/lel)]- 1 (P ion (t)) (solid 
line) and asymptotic behaviour according to Eq,(fl9|) 
(dashed) and Eq.(^) (long dashed). 

In the extreme opposite limit of large values of t we 
obtain the relations 



/(T) ^ r( I )r3/4 (T>>1) 



(47) and 

(48) 
(49) 



(pgS)) = 



n\) 



99 yu " tt Vr|/3 2 i 7 6 
r 8 

arctan 



4 
3^ 



2 r(i/4) 



1/4 



(50) 

(51) 
(52) 



/2 7 6i\3/4 



In this case the power law decays which characterize the 
diffusive dynamics of the excited Rydberg electron differ 
from the corresponding results of the PDM significantly. 
Even the characteristic exponents are changed. Accord- 
ing to Fig.|^ this dynamical regime is realized for interac- 
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tion times t which fulfill the relation i_i/4 <t <^t c with 
t_i /4 = T J p 2 /{2b 1 ). 

V. FULL MASTER-EQUATION 

Due to their simplicity and their applicability to all 
types of laser spectra the DCA rate equations are ideal 
for understanding the dynamics of Rydberg systems 
in the case of large laser bandwidths (compare with 
Eqs.(p^),([l9|) and However the DCA approxima- 

tion is not capable of describing coherent aspects of the 
laser-induced excitation process. In order to investigate 
the limits of applicability of the DCA rate equations in 
this section a more general approach is developed which 
is also capable of describing all coherent aspects of the 
excitation process. For the sake of simplicity we shall 
restrict our subsequent discussion to the case of phase 
fluctuations of the exciting laser field which deviate only 
slightly from a Lorcntzian spectrum and which can be 
modelled by Eqs.(|), (|5|) and (||). For these type of laser 
fluctuations we shall derive an approximate master equa- 
tion involving the density matrix elements of the excited 
Rydberg electron which are averaged over the fluctua- 
tions of the laser field. This procedure is a generalization 
of previous approaches which so far have been applied 
to atomic few level system only ||. In the special case 
of the PDM these subsequently derived density matrix 
equations reduce to our previous results of Ref. || . 

We start from the Schrodinger equation with Hamilto- 
nian ([lO]) with a fluctuating laser field as given by Eqs. 
(||) and (0). For the effective density operator 

P (t)= j2 \k)m\m)(m\3)e im)+wt){5ie - Ska) 



k,jGg,n 



we obtain the equation of motion 



P(t) = 

The self-adjoint Hamiltonian 
H dr . = \ n )( n \ e n + \g)(g\e - 



H dr ., P (t) -i<p(t) \g)(gl P (t) 



(53) 



(54) 



s J2(d ng \n)(g\ +h.c.) 



(55) 

describes the dynamics of the Rydberg system in the ab- 
sence of phase fluctuations and the stochastic process 
(j)(t) is defined by Eq.(||). In order to average Eq.(|5^) 
over all possible realizations of the stochastic process 0(7) 
it is convenient to introduce the averaged operators 



»(t) - 



0\ n/2 (-»)' 



d4Q n {<t>)p{t)p{<t>,t) (56) 



with n = 0, 1,2,3, ... . The (conditional) probability dis- 
tribution p(cj), t) obeys the Fokker-Planck equation 



[-+C]p{<t>,t) = 

with the Fokker-Planck operator 
d d 2 



(57) 



(58) 



This equation has to be solved with the initial condition 



p(<j>,0) 



V2^F 6XP V 2/36 



(59) 



which represents the stationary solution of the Fokker- 
Planck equation. The quantities 



)»W = H n (^[2 n+1 n\0b] 1/2 ) 



(60) 



with the Hermite polynomials H n are eigenfunctions of 
the adjoined Fokker-Planck operator C' with eigenvalues 
A n = n/3 JL3[ . Starting from Eq.(54) we obtain a set of 
coupled differential equations for the operators p( n '(t), 
namely 



H dr .,p (n \t) -n0pW(t) 



(61) 



+ b(n + l) [| 5 )<5l, (*)] - £ [\g)(g\, P (n - 1] (t) 

These equations have to solved subject to the initial con- 
dition 



pW(o) = <W(o). 



(62) 



According to Eqs.(|56]) and (60) p^(t) is the required 
density operator which is averaged over the phase fluctu- 
ations of the laser field. 

We may derive an approximate equation of motion for 
the averaged density operator p^'(t). As a first step we 
Laplace transform Eqs. (|6l|) , i.e. 



(63) 



-i\H dr .,pW(z)] + 



\g)(g\ 1 b{n+l)p {n+1) {z)-0p^- 1 Hz) 



In view of the large values of we are interested in 
we may neglect terms containing £o i n comparison with 
terms containing 0. Thus with the definition 



5?(*) 



(64) 



we arrive at the recursion relations 



n + 1 + -(ej - e- z) 



(65) 
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From these recursion relations we find 

*0' - 1 



1 + i 



2b/0 



3 



. e k - e - z 



36//? 



(66) 

Using Eqs.(|6~4|) and ( |66[ ) we may now eliminate p^\z) in 
Eq.(p3[) for n — 0. Performing the Laplace back trans- 
formation ( p7| ) and using the definition (p(t)) — p^(t) 
we finally obtain the master equation 



(P (*)> 



H dr .,{p{t)) 



(67) 



dr 



with the memory function 

f oo+zO 



^{ a fe (r)|fc)(fc|(p(t-r)>|. 9 )( g |+h.c.} 



1 
2^ 



^e- lzr a^(z) 



(68) 



— oc+iO 



In the limit of the PDM, i.e. for (3 — ► oo, this master 
equation reduces to the well known form [fL2| 



(P (*)> = -» H dr .,(p(t)) 



(69) 



i{[i,{p(t))it] + [ i(p(t )),it]} 



with the Lindblad operator 

L = V2b\g){g\. 



(70) 



VI. NUMERICAL RESULTS 

In this section numerical solutions of the master equa- 
tion ( |67j ) are compared with the corresponding solutions 
of the DC A rate equations (fl4|),(p3|) and ( p 5 |) . Details of 
the numerical technique for solving Eq.(|6J) are summa- 
rized in the appendix. On the basis of these comparisons 
the validity conditions for the applicability of the DCA 
and its accuracy can be tested. For this purpose we con- 
sider the laser excitation of a Rydberg system which can 
be described by quantum defect theory in a one-channel 
approximation (compare with Eqs.(O)). Typically this 
is a good approximation for Alkali atoms. 

The time evolution of the mean initial state probability 
(p gg (t)) and of the mean ionization probability (Pi on (t)) 
are depicted in Fig.|| for excitation at and well below the 
ionization threshold for different values of (3. In both 
cases it is assumed that the exciting laser field has a well 
defined amplitude and a fluctuating phase. 

Let us first turn to the case depicted in Figj^a: The 
spectrum of the laser field is close to Lorentzian (/3 ^> b), 



so that the asymptotic form Eq. (^|) applies well. Thus the 
effective bandwidth B as defined by Eq.(|l6|) is approxi- 
mately equal to the parameter b which characterizes the 
spectrum of Eq.(|]) and the parameter f3 might be in- 
terpreted as an effective cut-off frequency of the laser 
spectrum. Rydberg states are excited by the fluctuat- 
ing laser field well below the ionization threshold. The 
mean excited energy corresponds to a quantum number 
n res = (-2?)- 1 / 2 = 200. The laser bandwidth b and the 
laser-induced rate 7 are so small that the excited Rydberg 
states are located well below threshold, i.e. — e 3> b, 7. 
However, the values of b and 7 are large enough so that 
more than one Rydberg state around energy e is affected 
significantly by the laser field, i.e. Tfy,T^b > 1. The 
three curves of Fig.||a (solid, dashed and long dashed) 
correspond to different values of the effective cut-off fre- 
quency (3 of the laser spectrum. As many excited states 
are involved in the depletion of state \g) the initial stage 
of the time evolution is governed by an approximate ex- 
ponential decay of state \g) with rate 7 0]. This initial 
stage of the time evolution is independent of the fluctu- 
ations of the laser field. At larger interaction times with 
t > I/7 a coherent oscillation starts to appear in (p gg {t)) 
with the classical Kepler period Tj. This oscillation re- 
flects the time evolution of the electronic Rydberg wave 
packet which has been prepared by the fast depletion of 
the initial state \g). With each return to the core region 
this Rydberg wave packet might undergo a transition to 
state \g) thus increasing (p gg (t)). These coherent oscil- 
lations cannot be described by the DCA rate equations. 
However, due to laser fluctuations after a few Kepler peri- 
ods these coherent oscillations are damped out and merge 
into diffusive dynamics which is characterized by power 
law decay of the initial state \g). From this time on the 
dynamics of the Rydberg system under the influence of 
the fluctuating laser field is well described by the DCA 
rate equations. This is apparent from Fig.^a by compar- 
ing the numerical solutions of the master equation (solid, 
dashed and long dashed curves) with the asymptotic solu- 
tions of the DCA rate equations (circles and thin dashed 
curves ). According to the discussion presented in Sec. 
IV B 2 this diffusive dynamical behaviour appears when 



all coherent effects are damped out and disappears again 
at interaction times t > t c at which stochastic ionization 
starts to dominate. Physically speaking for these inter- 
mediate interaction times the excited electronic Rydberg 
wave packet starts to diffuse in energy space towards the 
ionization threshold. It reaches the ionization threshold 
roughly at time t c at which the ionization probability 
rises significantly from vanishingly small values to values 
close to unity. The early stages of this diffusion towards 
the ionization threshold are governed by the power law 
decay of Eq.(ft8|) for (p gg (t)) which is characterized by 
the exponent (—1/2). In the case of laser fluctuations 
which can be described by the PDM to a good degree 
of approximation this power law decay governs the time 
evolution up to the stochastic ionization time t c . How- 
ever, for non-Lorentzian spectra this is no longer the 



9 



case. In cases in which the non-Lorentzian effects can 
be described by the spectrum of Eq.(^|) the character- 
istic exponent of this power law decay is changed to a 
value of (—1/4) as soon as the interaction times become 
larger than the characteristic time t-i/A — Ti(3 2 / '(267) 
provided t_iu < t c (compare with Eq.(|5l|)). This non- 
Lorentzian effect is clearly apparent in Fig.pk where the 
characteristic times f_i/4 are indicated for b/(3 — 0.03 
and b/(3 — 0.2. With increasing values of j3 this char- 
acteristic time increases and this cross over phenomenon 



disappears for sufficiently large values of j3 as soon as 
t_i/4 > t c . At interaction times exceeding the stochastic 
ionization time t c the excited Rydberg wave packet has 
already reached the ionization threshold and the mean 
ionization probability rises to a value of unity. This 
asymptotic long time behaviour of the excitation dynam- 
ics is described to a good degree of approximation by 
Eqs.([32j) and (^3|). That is apparent from a comparison 
of the thin dashed lines of Fig.||a with the corresponding 
numerical solutions of the master equation (|67|). 




FIG. 3. Mean initial state probability {p gg (t)) and mean ionization probability (Pion(t)) as a function of interaction 
time t (in units of I/7) for different values of /3. 

a) excitation well below threshold : 6/7 = 5, Try = 2,n res = (-2e)~ 1/2 = 200, a = 0.1. PDM-limit b/j3 = 0: (solid 
curve), b[B = 0.03: (dashed curve), b/f3 = 0.2: (long dashed curve), homogeneously spaced energy level limit according 
to Eqs.(^J,^): (circles) and long time estimates according to Eqs.(^,^): (thin long dashed curves) 

b) excitation at threshold : 6/7 = 120, Try = 10,e/ 7 = -63.27, (-2e)" 1/2 = 200, a = 0.1. PDM limit (b/(3 = 0): 
(solid curve), strongly non-Lorentzian situations b//3 = 3: (dashed curve) and b/(3 = 10: (long dashed curve). Long time 
estimates according to Eqs.(B2,B3): (thin long dashed curves). 
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In Fig. [pa the laser bandwidth is so large that the sig- 
nificantly excited energy region [e — B, e + B] (before the 
onset of the electronic diffusion process) contains already 
the ionization threshold. Thus the excited Rydberg sys- 
tem is ionized significantly already in the early stages 
of the time evolution. As B > 7 this early stage of 
the ionization process is well described by the DCA rate 
equations which yield (compare to Eqs.(p9|) and (|40|)) 



(P.. 



(P, 



ft 



gg(t)) = e 

on(t)) = -(1 

7 



-ft 



(71) 



These approximate solutions are obtained from Eqs.(|14|) 
and (|2^) by neglecting (p n n(t)) in comparison with 
(pgg{t)). According to Eqs.(p9|) and (40) this initial 



ionization process saturates as soon as the mean ini- 
tial state probability and the mean ionization probability 
have reached the values Ag p and r/7. At these interac- 
tion times we still have jAs p t <C 1. These character- 
istic aspects of the laser-induced excitation process are 
clearly apparent in Fig]j|b. Physically speaking in this 
initial stage of the excitation process the Rydberg elec- 
tron is ionized with probability r/7. With a probability 
of (1 — r/7) an excited electronic Rydberg wave packet 
is prepared after a time of the order of 1/7. This wave 
packet is formed by a coherent superposition of all Ryd- 
berg states within the dominantly excited energy inter- 
val [e — B, 0]. Depending on the actual value of the laser 
bandwidth the coherent dynamics of this electronic wave 
packet is damped sooner or later. After the destruction 
of all coherences to a good degree of approximation the 
subsequent dynamics is governed by the DCA rate equa- 
tions. In Fig. ||3 small coherence oscillations are visible 
in the time evolution of {p gg (t)). As soon as the interac- 
tion time exceeds the stochastic ionization time the ex- 
cited Rydberg electron starts to ionize significantly. The 
time evolution of this stochastic ionization process is well 
described by Eqs.(|32|) and ( p3| ) within the framework of 
the DCA rate equations. 



VII. SUMMARY AND CONCLUSION 



The dynamics of an electronic Rydberg wave packet 
under the influence of a fluctuating cw-laser field has 
been discussed. It has been shown that for large laser 
bandwidths its dynamics can be described by Pauli-type 
rate equations for the relevant density matrix elements 
of the excited Rydberg electron averaged over the laser 
fluctuations. These rate equations are valid for arbitrary 
types of laser fluctuations and their dynamics is deter- 
mined by the spectrum of the laser field only and not by 
any of the higher order correlation functions. The valid- 
ity of these rate equations has been investigated in detail 
for a special class of phase fluctuations of the laser field. 

With the help of these rate equations we have investi- 
gated the dynamics of a laser excited Rydberg electron 
for long interaction times. At these interaction times 
the dynamics of the Rydberg electron are dominated by 
stochastic diffusion in energy space towards the ioniza- 
tion threshold which leads finally to stochastic ionization. 
This diffusion process is accompanied by a characteris- 
tic scenario of power law decays. Analytical expressions 
have been derived for these power laws and their associ- 
ated characteristic exponents. These analytical expres- 
sions exhibit in a clear way that the asymptotic power 
laws are independent of the quantum defect of the ex- 
cited Rydberg states and to which extent they depend on 
details of the laser spectrum. In particular, it has been 
demonstrated that the characteristic exponents which de- 
scribe the process of stochastic ionization are completely 
independent of the laser spectrum. However, the initial 
stages of the diffusion of the excited Rydberg electron 
depend on details of the laser spectrum. 

Support by the Deutsche Forschungsgemeinschaft 
within the SPP 'Zeitabhangige Phanomene und Meth- 
oden' is acknowledged. 



APPENDIX A: NUMERICAL SOLUTION 
TECHNIQUE FOR THE MASTER EQUATION 

In this appendix an efficient numerical method for solv- 
ing the master equation ( |67j ) is outlined. 

First, we make a rearrangement of Eq.(|67j) splitting it 
into a PDM-part and an additional part that disappears 
in the limit (3 — ► 00, namely 



P(t) 



H effP (t)-p(t)Hl f} ] +L P (t). (Al) 



Thereby we have introduced the effective non-Hermitian 
Hamiltonian 



H eff = H dr .-ib\g)(g\, (A2) 
the damping operator 

Lp(t) = 2b\g)(g\p(t)\g)(g\ (A3) 
+ b f dTj2{\n)(n\p(t)\g)(g\w n (T) +h.c.} 

and a memory function 



w n = lim fle- nT - a° (r). 

£1 — *oo 



(A4) 
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In the PDM-lim it, i .e. — ► oo, the second term on the 
right side of Eq.(A3) disappears and th e m aster equation 
reduces to Eq.(|69|). Integration of Eq.(Al) yields 



P {t) = U(t) P (0)U\t) + (A5) 

) 

dt'e(t - t')U{t - t')Lp(t')Ul (t - t') 



where U(t) — exp[— iH e fft] is a non unitary time evo- 
lution operator. If the initial condition is taken to be 
p(0) = \g)(g\, the Laplace transformed matrix elements of 
the density operator pij(z) = C z pij(t) = / °° dte lzt pij (t) 
become 



+ 



plk(z) = Hl ggk (z) 1 + 2bp gg (z) 

b^2^Hi ngk (z)w n (z)p n g(z) + Hig nk (z)w* n (-z)p gn (z)Y 

H a bcd{z) 

1 

2^ 



(A6) 



(A7) 



d Zl C Zl+t0 (a\ U{t) \b) [C Z1 - Z+M (d\ U(t) \c) 



with w n (z) = 1 — <5° (z). The Laplace transform ed t ran- 
sition amplitudes C z (i\ U(t) \j) appearing in Eq.( A7) are 
easily calculated 0||. 



C z ( 9 \U(t)\g)= _ 

z + ib — e — 2->{z) 

C z (g\ U(t) |n) = C z (n\ U{t) \g) 

— iEgdng 



(A8) 
(A9) 



(z - e n )[z + ib—e — ' 



Cz (n\U(t)\l) = 
z-e n l 

where = Y.n+g 



(A10) 



(z — e{)[z + ib — e — £(2 



\egd„ 



is the self energy of state \g) 
and the Kronecker-symbol S n i turns into a Dirac-delta 
function 5(e n —ei) for energy normalized continuum states 
|e„) and |e/). For the subsequent treatment it is conve- 
nient to introduce the expressions 



a n (z) 

Pn(z) 

E n (z) 



Png(z) 



[1 + 2bp gg {z)\e d ng , 

Pgn(z) 
[1 + 2bp gg (z)]e d ng ' 

Hgngg(z) p , . _ Hgg n g 



£od 



rig 



£od 



(All) 
(A12) 
(A13) 



a g 



Thus Eq. (|Aq) and Laplace transformation of Eq.(Al) 
yield 



Pgg( z ) 



K(z) 



l-2WC(z)' 



(A14) 
(A15) 



1 



2n(z + iO)(l-2bK,(z)) J 
with the characteristic kernel 



de[a e (z) - (3 e (z)] 



K(z) = H gggg {z) + b(j2+ J dn{e]}\e Q d n g\ 2 {E n {z)w n {z)a n {z) + F n (z)w*(-z)p n (z)} (A16) 



and with 



ai(z) 



I3i{z) 



™;(-z)C^)Mz g | 2 (>^) + + (1 - Ji(z)) (e^z) + T,(z)) 

(i - g,(*))(i - Ji{z)) - (^(zKi-^FQW 2 !^! 4 

Wl [z)d (z)\e dig\ 2 (Ei {z)+T l {z)) + (l-G l {z)) (f, (z) + St (z) 
(1 - Gz(*))(l - JiW) - {Mz^i-zWdizY^ 



(A17) 



(A18) 



The non-diagonal couplings (ai, Pi «-> a n , (3 n , n ^ I) due to the sum in Eq. ( |A6| ) give rise to the expressions Si(z) 
and Ti(z) appearing in Eqs.(Al7) and (A18), namely 



Si(z) = b{ +J q dn(e)}\e d ng \ 2 [f3 n (z) 



F n {z)-F l {z) F l {z)-E n {z) 

- Oi n (z) 



n^{g,l} 



Q - e r , 



T l (z) = b{ ]T +J q dn(e)}\e d„ 



Ej^Eijz) 
a n (z) p n {z) 



e n - e; - z— 2i0 
E t (z)-F n (z) 



Q - £r 



The diagonal couplings of the a n and P n yield Ji{z),Gi{z) and Ci(z), i.e. 



e n - ei+ z + 2i0 



(A19) 
(A20) 



12 



z + 2iQ 

bwi(z) f°° dz 1 {l + e(-ei)\e di g \ 2 [z 1 +ib-I-T l (z 1 +iQ)]~ 1 (z 1 ~ei+iO)- 1 } 
<,!i ' ' ~ " - ' ~ " ib-e- - z - iO)] {z± - e t + iO) ' 



2tt y_oo [zi - 

bwf(-z) f°° dzi { 1 + e (- e i)\ £ odi g \ 2 [zi 

2n [zi+ib-e—E(zi+iO)](zi — z 



-ib 



T,(zi — z — iO)] {z\ — z ■ 
e ( - iO) 



-zO)- 1 } 



(A21) 
(A22) 

(A23) 



All information about the quantum system is con- 
taine d in the self energy In order to solve 

Eqs.( [A17l - [A20|) for a given value of z, the coefficients 



Hg g gg, E[, F[, J[ and G[ have to be calculated. Using the 
energy levels and dipole matrix elements of Eqs.(|ll]) the 
self energy becomes 



with 



S(z) = Slu — i^- + iy- ; — r — - 

v ' 2 ' 1 - exp(-2mv(z))] 

v{z) = {-2z)- 1 ' 2 + a 



(A24) 
(A25) 



and with the (non-resonant) quadratic Stark-shift contri- 
bution Suj. In Q we calculated the quantity Hgggg (f(z) 
in that work) by contour integration. In an analogous 
way the quantities Ei,Fi, Ji and Gi can be calculated but 
for sake of brevity we do not give them here explicitly. 
Starting with S n = T n = 0, Eqs.( |AT4[|A"20| ) are solved 
by iteration. Actually we found the non-diagonal cou- 
pling terms Si,Ti to be very small in comparison with 
the diagonal couplings so that this iteration converges 
very rapidly. 
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